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We point out that cosmological constraint from the axion isocnrvature perturbation is relaxed if 
the Higgs field obtains a large field value during inflation in the DFSZ axion model. This scenario is 
consistent with the Higgs inflation model, in which two Higgs doublets have non-minimal couplings 
and play a role of inflaton. 


I. INTRODUCTION 

Strong CP problem is one of the remaining myster¬ 
ies of the standard model (SM). Among solutions to the 
strong CP problem proposed so far, the axion solution [l| 
is the most attractive one 0. In axion models, the global 
U(l) symmetry, which is called Peccei-Quinn (PQ) sym¬ 
metry, is spontaneously broken and the resulting pseudo 
Nambu-Goldstone (NG) boson takes a roll of axion to 
solve the strong GP problem. However, cosmological ef¬ 
fects of the axion are highly non-trivial. If the PQ sym¬ 
metry is broken after inflation, topological defects are 
formed and they may be harmful for cosmology. On the 
other hand, if the PQ symmetry is broken during infla¬ 
tion, the axion obtains quantum fluctuations which result 
in the cold dark matter (GDM) isocurvature perturba¬ 
tion Q. 

In the standard cosmological scenario in which the PQ 
symmetry is already broken during inflation and the PQ 
scalar (j) is stabilized at the potential minimum \(j)\ = 
fa/ with fa being the PQ scale, the magnitude of the 
axion isocurvature perturbation is given by 

a 2Trfa9i ’ 

where Hint denotes the Hubble scale during inflation and 
0i is the initial misalignment angle. The CDM isocur¬ 
vature perturbation is given by ^cdm = ra(26a/a) with 
To being the axion fraction in the present CDM energy 
density, and is evaluated as 

^CDM - 1 X (Y^T4iv) (lO?Sv)- 

The Planck constraint on the uncorrelated CDM isocur¬ 
vature perturbation reads S'cdm ^ 1-4 x I0“® 0. For 
the axion window 10® GeV ^ fa ^ 10^® GeV, the infla¬ 
tion scale is constrained to be iJinf ^ GeV taking 

account of non-Gaussian fluctuations [^ @|. In particu¬ 
lar, if we demand that the axion is a dominant compo¬ 
nent of GDM, the constraint reads Hini ^ 10^ GeV. It 
excludes high-scale inflation models unless PQ symmetry 
is restored after inflation. But in a class of axion models 
with domain wall number larger than one, such as the 
DFSZ axion model [ 3 ], the formation of axionic domain 
wall leads to a cosmological disaster. 


There are some ways to solve or relax the axion isocur¬ 
vature problem. For example, the PQ scale during infla¬ 
tion may be much larger than that in the present vac¬ 
uum. In this case, the isocurvature perturbation is sup¬ 
pressed by the ratio of the PQ scale at present and dur¬ 
ing inflation Q : see also [9l-[l^ . Another way is to make 
the axion heavy during inflation so that it does not ob¬ 
tain large scale quantum fluctuations due to the stronger 
QGD (la - fisj l or the explicit PQ breaking term [I^ [l3 |. 
The non-minimal kinetic term of the axion may also relax 
the constraint [l^- 

In this letter we consider a variant type of the sce¬ 
nario that suppresses the axion isocurvature perturba¬ 
tion. To this end, we consider the DFSZ axion model in 
which there are two Higgs doublets with PQ charges. If 
the Higgs bosons have large field values during inflation, 
the effective PQ scale is much larger than the present 
vacuum; the angular component of PQ scalar is mas¬ 
sive during inflation and does not obtain quantum fluc¬ 
tuations while the massless mode mostly consists of the 
pseudo scalar Higgs, which later becomes massive and 
does not lead to observable isocurvature perturbations. 
Actually, the effect of small mixing between the pseudo 
scalar Higgs and PQ scalar leads to the axion isocurva¬ 
ture perturbation, which is suppressed by the large field 
value of Higgs during inflation.®' In particular, we will 
show that it can be consistent with the Higgs inflation 
scenario [l^ in which the Higgs boson takes a roll of 
inflaton due to the non-minimal coupling to gravity. 


II. TOY MODEL ANALYSIS 

Before going into the detailed study, we first roughly 
sketch the basic idea in a toy model. Let us consider the 
model with two gauge singlet complex scalars (j) and S, 
both have opposite PQ charges. Then we write down the 
Lagrangian as 

£ = - \d^s\^ -V- X{c/^S^ + h.c.), (3) 


^ A similar idea has been proposed very recently in Ref. [T3| . 
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where A can be taken real and positive without loss of 
generality, and 

V = X<^ -y) + + As| 5|4 + X^s\^f\Sf. 

(4) 

Here A^, Xs and rrig are real and positive, rj is defined 
as rj = A^Dw/a with TVdw being the domain wall number 
which is a model dependent integer. In the following, 
we consider the case A^s < 2A for simplicity. For the 
stability of the potential at large IS”! and \(j)\, we need 


Therefore, if this exceeds the Hubble scale during infla¬ 
tion, 00 does not develop quantum fluctuations. This 
condition is written as 

A>(^)’, ( 11 ) 

But 00 constitutes a small fraction v^/vs) of massless 
mode o, hence the fluctuation of O 0 during inflation is 
estimated as 


(l5a0)inf 


^0 Hjnf 

VS 27r 


( 12 ) 


Xci,Xs > (A — A0s/2)^. (5) 

A term proportional to (fiS can be forbidden by a Z 2 
symmetry, under which S transforms as S —>• —S. The 
vacuum lies at S' = 0 and \cj)\ = rjjV^ where the NG 
mode, axion, is exactly the angular component of (j>'. 


V / .00 

—pz exp i — 

^ V ^ 


( 6 ) 


However, the role of (j) and S are inverted if S has a 
large field value during inflation. Actually it can be as 
large as ns = |S| ^ iSinf/v^As and S can be the infia- 
ton consistent with Planck observations for As ^ 10“^^ 
by introducing higher order non-renormalizable terms to 
make the S potential polynomial ( 2 ^ or for As 0(1) 
by adding non-minimal coup ling s to gravity |19l| or non- 
minimal kinetic terms of S 2l|. Let us write (p and S 
as 



In the vacuum, = rj and vs = 0. But they need not 
coincide with the vacuum values in the early universe. 
We assume that initially S has a large field value vs- 
Then v^ is given by 


n 0 ~ max 


•n, 


l2X — A0S 
2Xrk 


vs 


( 8 ) 


If we introduce a Hubble induced mass term for (p ns 
V ~ with c being a constant, it can also sta¬ 

bilize the PQ scalar at around ^0 ^ cH/ ^/X^. The mass¬ 
less mode, which we denote by a, mostly consists of the 
angular component of S ii vs ^ v^: 

a OS as - 00 . (9) 

vs 

It is mostly the pseudo scalar as that obtains quantum 
fluctuations of ~ iLinf/27r. On the other hand, the or¬ 
thogonal component, mostly consisting of a^, obtains a 
mass of 


- Xvl- 


( 10 ) 


Thus we obtain the axion isocurvature fluctuation at the 
QCD phase transition as 

(^) == 3 %-^ m 

\ / QCD 

Compared with the standard scenario with the broken 
PQ symmetry during inflation, the isocurvature pertur¬ 
bation is suppressed by the factor ~ fa/ vs “C 1. Since 
we do not introduce an explicit PQ breaking, there is 
always a NG mode, but the massless mode during infla¬ 
tion needs not coincide with that in the present vacuum. 
The axion at present was massive during inflation, while 
the massive mode was massless. The CDM isocurvature 
perturbation then is given by 


5eo„.lxl0-..(^)““(^)(A). 

(14) 

Thus inflation scale as large as Hmi ^ 10^^ GeV may be 
allowed for vs ~ 10^® GeV even if the axion is a dominant 
component of CDM. 

In this class of scenario, S has a large field value ini¬ 
tially and hence it begins a coherent oscillation after in¬ 
flation and must decay into SM particles. Later we will 
see that S itself can be the infiaton. The oscillation of 
S potentially causes resonant particle production of ax- 
ions [ 2 ^ and may induce the dynamical motion of (j)- We 
denote the frequency of S after the onset of the oscilla¬ 
tion as First, we consider the axion production due 

to the oscillation of S. If the condition 


VXvs 


< 1 , 


(15) 


holds, the axion production of the broad resonance type 
does not occur even just after the onset of the oscillation 
of S. In such a case, the narrow resonance would harm¬ 
fully produce axion particles. However, the narrow reso¬ 
nance is ineffective if the decay rate of S is sizable [^ . 
Second, let us consider the motion of (j) caused by the 
oscillation of S. The potential of (j) varies as S oscillates 
and (j) may also oscillate around the origin. In such a case, 
the nonthermal restoration of the PQ symmetry [23 - [T7I| 
could occur due to the parametric resonance at the ori¬ 
gin because of the self interaction of (j). The parametric 
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resonance at the origin becomes most effective when S 
decays just after the onset of oscillation of S. In such a 
case, the condition for the PQ symmetry not to be re¬ 
stored is obtained from the condition that the number of 
(j) oscillation within one Hubble time remains less than 
^ 0(10^) until the oscillation amplitude of </> reduces to 
~ fa [13 ■ This leads to a constraint 

n 

<0(10^). (16) 

where n = 1/2 (1) for the matter (radiation) dominated 
background evolution, where we have assumed the in¬ 
equality dm is marginally satisfied: \/\vs ~ Hmt- 
Next, let us consider the thermal effects on fi. In order 
for the PQ symmetry not to be restored after inflation 
due to thermal effects, we need 

Xi,s<xJ7ff^) , ( 17 ) 

\ max / 

where T^ax denotes the maximum temperature after in¬ 
flation and assumed that S is in thermal equilibrium due 
to interactions with SM particles. Otherwise, the PQ 
symmetry may be restored after inflation and isocurva¬ 
ture constraint disappears, but there arises a serious do¬ 
main wall problem. Conditions (fTTl) . (ITSl) . (fTO and m 
are satisfied for reasonable choice of parameters. 



III. DFSZ AXION MODEL 


scalar. Quantities with subscript J indicate those in the 
Jordan frame. Here is assumed to couple with up- 
type quarks and Hd with down-type quarks and charged 
leptons. This can be done by assigning the PQ charge, 
e.g., 1 to the PQ scalar /), —2 to and 2 to right- 
handed up-type quarks. In this model, the domain wall 
number TVdw is equal to 6. In the present vacuum, the 
PQ scalar obtains a VEV of |(()| = ry /^/2 and Higgs bosons 
have electroweak scale VEVs. Note that since the Higgs 
bosons also have PQ charges, the axion is a mixture of 
the angular components of (j) and Higgs bosons: 


CL — CL(p 


sin 2/3 
1 


{vdau + Vuttd) , 


where we have defined 


(23) 


rrO 

H„ = ^exp 




and tan/3 = Vujvd- However, the mixing angle is sup¬ 
pressed by the ratio Vu,d/fa and hence the massless mode 
almost consists of a^. Actually it is this small mixing 
that admits the axion-gluon-gluon coupling, required for 
solving the strong CP problem. Another massless mode, 
ac = cos Pad — sin /3a„, is eaten by the Z-boson and the 
other orthogonal combination becomes massive, which 
we identify as the pseudo scalar Higgs and denote by 
Qh- For large enough Vu,Vd ^ v^, however, the axion is 
identified as 

v± 

a ~ cos Pau + sin Pad -(25) 

vs\n2p 


Let us move to the DFSZ axion model, two Higgs dou¬ 
blets are introduced [l^. In this model, S'^ in the toy 
model may be identified with the gauge invariant 
combination of the Higgs doublets HuHd- The action 
in the Jordan frame is given by 

S = J (Cg +Lj-Vj + , (18) 


with V = + Actually, in the Higgs inflation 

model, Vu and Vd can take so large values that the axion 
isocurvature perturbation is suppressed as shown in the 
toy model in the previous section. 

The detailed analysis of Higg s inflation with two Higgs 
doublets are found in Ref. [^. Thanks to the PQ sym¬ 
metry, allowed terms are limited. After the conformal 
transformation where 


where 


^+UHu?+UHdl^Rj, 

(19) 

- \V^Hdl - 

(20) 


Vj = ml\Hu\‘^ -t ml\Hd\'^ + {Xp^HuHd + h.c.) 

+ Xu\H^\^ + Xd\Hdf 
+ X^d\Hu\'^\Hd\^ + X'^^\H^Hd\^ 

+ Xu4,\Huf\pf + Xd^\Hdf\pf + Vi\P\), ( 21 ) 

vm) = x^(^\p\^-'pj , ( 22 ) 

and contains kinetic terms for SM quarks, lep¬ 

tons and gauge bosons and yukawa couplings between 
quarks/leptons and Higgs bosons and Rj is the Ricci 


r>2 _ 1 I 

+ Ml + Ml ’ ^ 

the Einstein frame action becomes 

S = j + Ce-Ve+ > (27) 

with subscript E indicating the Einstein frame. The 
scalar potential is given by 

VEiH^,Hd,P) = ^. (28) 

Focusing on only terms with fourth powers of and hd, 
it is expressed as 


Ve XuV^ + Xdvj -I- Xudvlvj 
Ml~ ^{l + fuvl+pdvlY 
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where \ud = ^ud + Thus it is easy to see that 

the potential becomes flat for ^ Mp/y^ or Vd ^ 
Mp/\/^. Terms involving </) do not affect the inflaton 
dynamics as we will see. There is a stable inflationary 
path along 

'^u _ ‘^^d^u ~ ^ud^d ("301 

2A„^d — Xud^u ’ 

if ‘^Xd^u Xud^d ^ 0 and 2A^^d X^d^u ^ d- Then the 
potential energy for the inflaton is given by 

_ _ XuXd — X1^/A _ / ^ 2x/V6Mp\^ 

K 4(A,e^ + Ad^^ _ XudUd) V / ’ 

(31) 

where y is the canonically normalized held in the large 
held limit: 


X = 


Mp log 




(32) 


Thus it reduces to the single field Higgs inflation model. 
Assuming Xu ^ XdXud ~ A^^ and ^d, we need 


e - 5 X (33) 


to reproduce the density perturbation observed by 
Planck 0, where ^ and Aeff denote typical values of ^’s 
and A’s, respectively. To explain the 125 GeV Higgs bo¬ 
son, Aeff should be 0 ( 1 ), but the running tends to make 
A smaller at high energ y sc ale and hence ^ may be able 
to take a smaller value . 

Next, let us see the behavior of the PQ scalar during 
inflation. It obtains a VEV of 


v,p ~ max 


1, 



(34) 


where f/i(~ ~ t’d ~ 10 Mp/i/?) denotes the typical 

held value of the Higgs boson during inflation. The dy¬ 
namics after inflation is effectively the same as that stud¬ 
ied in the previous section. Note that obtains a mass 
of 


m 


2 



AM|, 


Thus we impose the condition 


A> 




inf 


M2 


(35) 


(36) 


to suppress the fluctuation of during inflation. Also 
note that the reheating temperature after inflation may 
be as high as ~ 10^^ GeV taking account of efficient par¬ 
ticle production and thermal dissipative effects |3a - K^ . 
To avoid the PQ symmetry restoration, we need a con¬ 
dition similar to (IT71) : 


Xucj }; Xd 


<Xd 


fa 

Tma, 


(37) 


The pseudo scalar Higgs, Uh, is almost massless during 
inflation and develops quantum fluctuations. Note that 
the canonically normarized held is a-h = auMp/\f^Vh 
which obtains a fluctuation of 5ah — i7inf/27r, while 
the mixing between and ah is given by ^ v^jvh- 
Then we find the effective fluctuation of as 5a^/'V(p ^ 
i/jMinf/(27rMp). Thus the axion isocurvature fluctuation 
is given by 


( ^ VfHiui 
\ / Qco ‘2'jTMpOi 


(38) 


The CDM isocurvature perturbation is given by (|T4)) with 
Vs Mp/y/^. Thus the constraint is significantly re¬ 
laxed. 


IV. DISCUSSION 

In this letter we have shown a way to suppress the ax¬ 
ion isocurvature perturbation. Even if the PQ scalar it¬ 
self does not have a field value much larger than fa during 
inflation, other PQ charged scalar fields can have large 
field values and it results in suppression of the isocurva¬ 
ture perturbation. Interestingly, the electroweak Higgs 
bosons may take such roles in the DFSZ axion model, 
and the Higgs bosons can also be the inflaton. 

So far, we have considered the quartic interaction term 
between the PQ scalar and the Higgs bosons XcjPHuHd + 
h.c.. There are other possible choices. For example, let 
us consider the case with the cubic term V = ^4>HuHd + 
h.c.. In this case, (j) settles around (j) ~ M/Aii<^,d 0 and has 
the mass of ~ Xu4,,d4>Vh during the inflation. After the 
inflation, cj) starts to oscillate. If p, < /a holds, S will not 
pass through the origin and the restoration of the PQ 
symmetry due to the particle production does not occur. 
The higher dimensional terms may also play important 
roles such as stabilization of the PQ scalar at some held 
value. 

Although we considered the case of Higgs boson as the 
inflaton, we can also regard the PQ scalar as the inflaton 
in a similar way, by introducing a non-minimal coupling 
as [HI'" 


C = ^^\ct>fRj. (39) 

The PQ scalar can take a role of inflaton for ^ 
5 X lO'^y/X^. By noting that the canonically normal¬ 
ized axion field at > Mp/yJ^ is given by (itj, = 
a^Mp/ and it obtains quantum fluctuations of 

— Minf/27r, we would have the axion isocurvature per¬ 
turbation of the same expression as (l38l) . if the PQ sym¬ 
metry is never restored thereafter. There is a subtlety. 


The PQ scalar as the inflaton in the context of running kinetic 
inflation was mentioned in Ref. Ea. 
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however, in this scenario. Because the canonically nor¬ 
malized PQ scalar has an oscillation amplitude of ^ Mp 
just after inflation, it induces efficient axion production 
when it passes through the origin of the potential, and 
the PQ symmetry may be restored nonthermally, lead¬ 
ing to formation of axion domain walls if the domain 
wall number is larger than one . We can arrange the 
model so that the Higgs bosons also have large VEVs 
during and after inflation just as we have done in the 
previous section for (j). Then the PQ symmetry is broken 
by Higgs VEVs even though (f> is nonthermally trapped 
around the origin. Even in such a case, however, the 
model has a Z 2 symmetry under which cj) ^ —cj) and any¬ 
thing else uncharged, hence the domain wall formation is 
not avoided when </> relaxes to the minimum. Again, in¬ 


troducing V = ^(j)HuHd h.c. in the potential, instead 
of X(j)'^HuHd + h.c., by assigning appropriate PQ charges, 
may help the situation but a complete analysis is beyond 
the scope of this paper. 


Acknowledgments 

This work was supported by Grant-in-Aid for Scien¬ 
tific research 26104009 (KN), 26247042 (KN), 26800121 
(KN), MEXT, Japan. The work of M.T. are supported 
in part by JSPS Research Fellowships for Young Scien¬ 
tists and by the Program for Leading Graduate Schools, 
MEXT, Japan. 


[1] R. D. Peccei and H. R. Quinn, Phys. Rev. Lett. 38, 1440 
(1977). 

[2] For a review, see J. E. Kim, Phys. Rept. 150, 1 (1987); 

J. E. Kim and G. Carosi, Rev. Mod. Phys. 82, 557 (2010) 
arXiv:0807.3125 [hep-ph]]. 

[3] For a review on axion cosmology, see M. Kawasaki and 

K. Nakayama, Ann. Rev. Nucl. Part. Sci. 63, 69 (2013) 
arXiv:1301.1123 [hep-ph]]. 

]4] P. A. R. Ade et al. [Planck Collaboration], 
arXiv:1502.02114 [astro-ph.CO]. 

[5] M. Kawasaki, K. Nakayama, T. Sekiguchi, T. Suyama 

and F. Takahashi, JCAP 0811, 019 (2008) 

arXiv:0808.0009 ]astro-ph]]; C. Hikage, M. Kawasaki, 
T. Sekiguchi and T. Takahashi, JCAP 1307, 007 (2013) 
arXiv:1211.1095 [astro-ph.CO]]. 

[6] T. Kobayashi, R. Kurematsu and F. Takahashi, JCAP 
1309, 032 (2013) ;arXiv:1304.0922 [hep-ph]]. 

[7] M. Dine, W. Fischler and M. Srednicki, Phys. Lett. B 
104, 199 (1981); A. R. Zhitnitsky, Sov. J. Nucl. Phys. 
31, 260 (1980) [Yad. Fiz. 31, 497 (1980)]. 

[8] A. D. Linde and D. H. Lyth, Phys. Lett. B 246, 353 
(1990); A. D. Linde, Phys. Lett. B 259, 38 (1991). 

[9] K. Nakayama and N. Yokozaki, JHEP 1211, 158 (2012) 
arXiv: 1204.5420 [hep-ph]]. 

[10] T. Moroi, K. Mukaida, K. Nakayama and M. Takimoto, 
JHEP 1411, 151 (2014) arXiv: 1407.7465 [hep-ph]]. 

[11] M. Fairbairn, R. Hogan and D. J. E. Marsh, Phys. Rev. 
D 91, no. 2, 023509 (2015) arXiv:1410.1752 [hep-ph]]. 

[12] K. Choi, E. J. Chun, S.' H. Im and K. S. Jeong, 
arXiv:1505.00306 [hep-ph]. 

[13] G. R. Dvali, hep-ph/9505253 

[14] K. Choi, H. B. Kim and J. E. Kim, Nucl. Phys. B 490, 
349 (1997) hep-ph/9606372 . 

[15] K. S. Jeong and F. Takahashi, Phys. Lett. B 727, 448 
(2013) arXiv: 1304.8131 [hep-ph]]. 

[16] M. Dine and A. Anisimov, JCAP 0507, 009 (2005) 
hep-ph/0405256 . 

[17] T. Higaki, K. S. Jeong and F. Takahashi, Phys. Lett. B 
734, 21 (2014) arXiv:1403.4186 [hep-ph]]. 

[18] S. Folkerts, C. Germani and J. Redondo, Phys. Lett. B 
728, 532 (2014) arXiv:1304.7270 [hep-ph]]. 

[19] F. L. Bezrukov and M. Shaposhnikov, Phys. Lett. B 659, 
703 (2008) arXiv:0710.3755 [hep-th]]. 


[20] K. Nakayama, F. Takahashi and T. T. Yanagida, 

Phys. Lett. B 725, 111 (2013) arXiv: 1303.7315 

[hep-ph]]; JCAP 1308, 038 (2013) arXiv: 1305.5099 
arXiv: 1305.5099 [hep-ph]]. 

[21] F. Takahashi, Phys. Lett. B 693, 140 (2010) 

arXiv:1006.2801 [hep-ph]]; K. Nakayama and F. Taka¬ 
hashi, JCAP 1011, 009 (2010) arXiv: 1008.2956 [hep- 
ph]]. 

[22] L. Kofman, A. D. Linde and A. A. Starobinsky, Phys. 
Rev. D 56, 3258 (1997) hep-ph/9704452 . 

[23] Y. Ema, R. Jinno, K. Mukaida and K. Nakayama, 
arXiv: 1504.07119 [gr-qc]. 

[24] L. Kofman, A. D. Linde and A. A. Starobinsky, Phys. 
Rev. Lett. 76, 1011 (1996) hep-th/9510119 . 

[25] S. Kasuya, M. Kawasaki and T. Yanagida, Phys. Lett. B 
409, 94 (1997) hep-ph/9608405 ; Phys. Lett. B 415, 117 
(1997) hep-ph/9709202 . 

[26] T. Moroi, K. Mukaida, K. Nakayama and M. Takimoto, 
JHEP 1306, 040 (2013) arXiv: 1304.6597 [hep-ph]]. 

[27] M. Kawasaki, T. T. Yanagida and K. Yoshino, JCAP 
1311, 030 (2013) arXiv:1305.5338 [hep-ph]]. 

[28] For a review on two Higgs doublet models, see 
G. C. Branco, P. M. Ferreira, L. Lavoura, M. N. Re- 
belo, M. Sher and J. P. Silva, Phys. Rept. 516, 1 (2012) 
arXiv:1106.0034 [hep-ph]]. 

[29] J. O. Gong, H. M. Lee and S. K. Kang, JHEP 1204, 128 
(2012) arXiv: 1202.0288 [hep-ph]]. 

[30] A. Salvio, Phys. Lett. B 727, 234 (2013) arXiv: 1308.2244 
[hep-ph]]. 

[31] Y. Hamada, H. Kawai, K. y. Oda and S. C. Park, Phys. 
Rev. D 91, no. 5, 053008 (2015) arXiv: 1408.4864 [hep- 
ph]]. 

[32] F. Bezrukov, D. Gorbunov and M. Shaposhnikov, JCAP 
0906, 029 (2009) arXiv:0812.3622 [hep-ph]]. 

[33] J. Garcia-Bellido, D. G. Figueroa and J. Rubio, Phys. 
Rev. D 79, 063531 (2009) arXiv:0812.4624 [hep-ph]]. 

[34] K. Mukaida and K. Nakayama, JCAP 1301, 017 

(2013) arXiv: 1208.3399 [hep-ph]]; JCAP 1303, 

002 (2013) arXiv:1212.4985 [hep-ph]]; K. Mukaida, 

K. Nakayama and M. Takimoto, JHEP 1312, 053 (2013) 
arXiv:1308.4394 [hep-ph]]. 



